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Abstract 



A light-front treatment for spherical nuclei is developed from a relativistic 
effective Lagrangian and employing the mean field approximation. Minimizing 
the nuclear minus momentum subject to the constraint that, in the rest frame, 
the expectation values of the plus- and minus-momenta are the same leads to 
a formalism in which rotational invariance is recovered. 
Light cone variables have long been used to analyze high energy experiments involving 

1 



nuclear and nucleon targets. The most well-known use appears in the parton model. The 
Bjorken variable x is the ratio x = k + /p + where k + = k° + k 3 is the plus- momentum of 
the struck quark and p + is the plus- momentum of the target. Quark distributions represent 
the probability that a quark has a plus-momentum fraction x. Our focus here is on the 
distribution functions which describe the plus-momentum carried by the nucleons within 
the nucleus. Such distributions, which depend on p + of the struck nucleon, are needed to 
analyze deep inelastic lepton- nucleus scattering Q, and also enter in the analysis of high- 
momentum transfer nuclear quasi-elastic reactions such as (e, e'), (e, e'p), (p, 2p) 0- ||. One 
may also obtain similar distributions for the nuclear mesons, and these enter in the analyses 
of nuclear Drell-Yan experiments ||. If the light front formalism is used, these distributions 
are simply related to the square of the nuclear ground state wave functions ||. Many 
nuclear high momentum transfer experiments are planned, so that it is necessary to derive 
a relativistic formulation which uses p + variables closely related to experiments, and which 
incorporates the full knowledge of nuclear dynamics. 

If the relevant nuclear wave functions depend on p + , the canonical spatial variable is 
x~ = x° — x 3 . This leaves x + = x° + x 3 to be used as a time variable, with the light 
front Hamiltonian (x + development operator) as P~ = P° — P 3 . These are the light front 
variables of Dirac 0. There are clear advantages in using these variables, but a principal 
problem arises because the use of x + as "time" and x~,x± as "space" involves the loss of 
manifest rotational invariance. This is especially important in nuclear physics because the 
understanding of magic numbers rests on the 2j + 1 degeneracy of single particle orbitals of 
good angular momentum. Our purpose is to show how rotational invariance is recovered in 
nuclear mean field theory. Once this is accomplished, light front dynamics can be a useful 
tool in a variety of nuclear physics calculations. 

We start by assuming a field-theoretic Lagrangian in which the nuclear constituents are 
nucleons ip scalar mesons <fr and vector mesons V^, = d^V" — d u V )1 . The Lagrangian is 
a standard one j5|, and we use the meson- nucleon coupling constants and meson masses of 
Horowitz and Serot || . The use of light front dynamics allows us to compute the necessary 
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nuclear meson distribution functions using the same formalism that yields the nuclear bind- 
ing energy and density. Previously, one of us constructed a light front mean-field treatment 
of nuclear matter [TIJ using the same Lagrangian. Vector mesons were found to carry 35% of 
the total plus- momentum (mass), with each vector meson carrying an infinitesimally small 
amount of plus-momentum. Does such a result hold for finite nuclei? 

Light front quantization of a given Lagrangian is a standard procedure || . One first uses 
the Lagrangian to derive field equations which enable the identification of dependent degrees 
and their elimination in favor of independent ones. The procedure for the present Lagrangian 



is described elsewhere |[L0| , so we discuss only the most relevant features concerning the 
nucleons. Two of the four degrees of freedom of this spin 1/2 particle entering in the 
usual Dirac spinor must be redundant. These are separated using the projection operators 



A± = j7 7 , with ip± = A±ip. The independent fermion degree of freedom is chosen as ip + . 
In general the equation for depends on the interaction V + in a complicated manner. This 
dependence is eliminated by using the transformation of Refs. ||11||. The net result is that 



a transformed vector meson field = — enters in the nucleonic field equation, 

while enters in the equation for the vector meson field. 

The energy momentum tensor is T^, and the volume integrals of components T H and 
T ++ are the Hamiltonian P~ and momentum P + . We find 

T+- = (V ± 0) • ( V ± 0) + m 2 ^ 2 - (V ± V + ) ■ (V ± V~) - m 2 v V + V- + 2^ (id~ - g v V~) 
T++ = d + (f)d + (f) + (V ± V + ) ■ (V ± V + ) + m 2 v V + V + + 2^ ] + id + ^ + . (1) 

We choose to work in the nuclear rest frame so that P± | 1 J r ) = and P + \ ty) = Ma \ ^/). 
In this frame, the light-front Schrodinger equation for the ground-state wave function | ty) 
reads 

M 2 -l- P 2 

P~ I *> = p + ~ I *> = M A | (2) 

In light-front dynamics the goal is to solve the equation (P + P~ — Pj_) | ty) = M\ \ 
with fixed values for Pj_ and P + , using a Fock-space basis. This is not possible for systems 



containing many nucleons, so we need to construct a light front version of the Hartree- 
Fock procedure in which the energy of a Slater determinant is minimized. Minimizing the 
expectation value of P~ doesn't work because one gets a zero of P for an infinite value 
of P + . Instead we minimize the expectation value of P~ subject to the condition that 
the expectation value of P + is equal to the expectation value of P~ . This is the same as 
minimizing the average of P~ and P + . To this end we define a light front Hamiltonian: 
Hlf = \ (P + + P~)- This is not the usual Hamiltonian because light-front quantization is 
used to define all of the operators that enter. 

The wave function | consists of a Slater determinant of nucleon fields | $) times a 
mesonic portion, and the mean field approximation is characterized by the replacements: 
cf) -> | | ^), -> | V» | We quantize the nucleon fields using 

^(x)=J2(x-,x ± \n)e-^ x+ / 2 b n , (3) 

n 

in which the variable x represents both the spatial variables x_l,x~~ and the spin, isospin 
indices, and b n obeys the usual anti-commutation relations. The meson fields are also func- 
tions of x~,x±. The single particle states | n) are defined by minimizing the nucleonic 
contribution to the expectation value of | H LF | \1>). The Slater determinant | $) is 
defined by allowing A nucleon states to be occupied. The results of the minimization are 

p- | n)+ = (id + + 2g v V~) \ n)+ + (a ± ■ (p ± - g v V ± ) + (3{M + g s 4>)) | n>_ , 
id + | ra)_ = (a ± ■ (p ± - g v V ± ) + f3(M + grfj) \ n)+ , (4) 

where | n)± — A± | n). Minimization of | H LF | \1>) with respect to the meson field 
variables leads to the meson field equations. We find 

. X 

dx 



-Vi - (-2^- 1 +mi\{*\ <f>(x) | M>> = -g sPs (x), (5 ) 



-Vi - (2^) + j (tt I ^(r) I *> = ^(s), (6) 
where p s (x) = | ^(x)ip(x) | and ^(x) = | %l)(x)^ip(x) \ ^/). 



We next discuss how Eq. (£|) is solved. The first step is to use a representation [12| in 
which | n)± are each represented as Pauli spinors. The light front quantization respects 
manifest rotational invariance for rotations about the z-axis. Thus each single-particle state 
has a good J z . We use a momentum representation for the longitudinal variable in which 
the values of p + take on those of a discrete set: p+ = (2m + l)vr/ (2L) where m > and L 
is a quantization length. This means that the nucleon wave functions have support only for 
p + > 0. This spectrum condition is a requirement for exact solutions for any theory. The 
coordinate space representation is used for the _L variables. Then we have 

U^(x ± )e i ^- 1 / 2 ^ 
L^(*±K (Jz+1/2)0 

in which the upper (lower) entrees of the Pauli spinor correspond to m s = l/2(— 1/2), and 
x± = (x±, 0). When Eq. (|7|) is used in Eq. (|), one finds that the equations do not depend on 
the magnitude of J z \ solutions for ± a given magnitude of J z are degenerate. The functions 



(Pm, x ± I n). 



(7) 



U m (x±) and L m (x±) are expanded in a basis of B-splines of degree five in x± JT3|. One starts 
with a guess for the scalar and vector potentials. Then one solves Eq. (|j). The solutions 
are used to generate new scalar and vector nucleon densities, which in turn generate new 
potentials via Eqs.(|5|) and (0). This iteration procedure is completed until the solutions are 
stable. 

If these solutions are to have any relevance at all, they should respect rotational invari- 
ance. The success in achieving this is examined in Table I, which gives our results for the 
spectrum of 40 Ca. We see that the J z = ±1/2 spectrum contains the eigenvalues of all 
states, since all states must have a J z = ±1/2 component. Furthermore, the essential fea- 
ture that the expected degeneracies among states with different values of J z are reproduced 
numerically. 

The results shown in Table 1 are obtained using a basis of 20 splines, a box size of 
2L = 24 fm, and 24 Fourier components in the expansion of the wavefunction. This value 
of L is large enough so that our results do not depend on it, and the number of terms in the 
expression for the density is enough to ensure that the densities are spherically symmetric 



TA\. Another feature is that the spectrum with p + > has no negative energy states, so 
that in using the LF method one is working in a basis of positive energy states only 

The single particle energies e n reported in Table 1 are given by e n — p~/2. The values of 
e„ are essentially the same as those of the ET formalism to within the expected numerical 
accuracy of our program, as are the nuclear densities. This equality is not mandated by 
spherical symmetry alone because the solutions in the equal time framework have non- 
vanishing components with negative values of p + . We may understand the near equality of 
single particle eigenvalues using an analytic argument. First we use a representation in which 
appears in the equation for the nucleon field. That is, define (x \ n)' = e~ t9vA ^(x \ n), 
with d + A = V°. Then multiply the equation for | n)' + by 7 + and the equation for | n)'_ by 
7~~. Adding the resulting two equations gives 

( 7 °(p; - 2 7 °g v V° - 7 3 (2p + - p-/2)) | n)' = 2( 7± ■ p± + M + g s ct>) \ n)' . (8) 

Solving Eq. (§) in coordinate space is expected to lead to solutions in which the spectrum con- 
dition is not respected exactly. Equation (|8|) may be converted into a manifestly rotationally 
invariant equation using the relation: x~ = —2z, so that p + = id + = 2i-£p — > — i-^. The 
operator p + acts as a p 3 operator, and the result (|8|) looks like the Dirac equation of the equal 
time (ET) formulation (of eigenvalue Pn/2) except for an offending term —p~/2 multiplying 
the 7 3 . This term may be eliminated by including a phase factor: (x \ n)' = e iPnZ ^ 2 (x \ u)et- 
The result is that (x \ n) et satisfies the standard Dirac equation of the equal time formu- 
lation. This means that the eigenvalues e n must be approximately the same in the two 
formulations. The net result of these transformations is that if one were to include p + < 
in the LF calculation then the eigenvalues e n in the LF and ET calculations would be 
identical and the LF wave function would be expressed as 

I r oo 

(p + ,x ± \n} = ^ dze- i{p+ - p - /2)z e i9 * A{x ^ z) (z J x ± \n) ET . (9) 
V 2n J-od 

The relation (|9|) tells us that (p + ,x± \ n) peaks at p + ~ M — gvV°, with a width of the 
order of the inverse of the radius of the entire nucleus. Therefore, neglecting p + < is only 



a minor approximation, which is the reason why our above LF calculation, with p + < 
suppressed, leads to eigenvalues e n that are approximately equal to the ET results. 

One can also read off from Eq. (0) that the influence of the vector potential is to remove 
plus-momentum from the nucleons. Furthermore, the large value of the nuclear radius causes 
the region of support to be very narrow, so that {p + , x± \ n) is very small for negative values 
of p + . 

The nuclear energy is obtained by adding the contributions to the expectation value of 
P~ + P + from the nucleons, and scalar and vector mesons. Then straightforward calculation 
leads to the result 

M A = - ~ / d 3 r [g s <j>{r)p s (r) + g v V°(r)p B (r)] , (10) 

n<F 

where r = (x±, z = —x~/2). The above expression is the result obtained in standard equal 
time calculations, and the nuclear binding energy we obtain is essentially the same as that 
of Ref. §. 

The previous results mean that we have reproduced the standard good results for nuclear 
binding energy and density using a formalism that allows us to compute nucleon and meson 
plus-momentum distributions used to analyze various high energy nuclear experiments. The 
probability that a nucleon Jn{p + ) has a momentum p + is given simply by 

MP + )=E /^i|(xi,p + |n) | 2 . (11) 

n<F J 

It is convenient to express the distribution in terms of the dimensionless variable y = 
P + /(Ma/A). The result is shown in Fig. 1. The peak of the distribution is at y « 0.76, 
corresponding to p + of 710 MeV (which is the plus-momentum per nucleon). Our numer- 
ical calculations show that the average value of p + is 672.56 MeV (< y >= 0.72) out of a 
total expectation value P + /A of 930.56 MeV. (P + /A is the total energy per nucleon.) The 
distribution is not symmetric about its average value, as it would be if a simple Fermi gas 
model were used. Both of these effects are caused by the presence of nuclear mesons, which 
carry the remainder of the plus-momentum. 



The meson distributions are computed using 

f s {k + ) = J d 2 k ± | (* | o(k) | *> | 2 (12) 
for the scalar meson distribution, and 

f v {k + ) = fd 2 k ± £ | | a(k,u) | | 2 (13) 

w=l,3 

for the vector meson distribution. The operator a(k) destroys a scalar meson of momentum 
k = (k + , k±), and a(k, destroys a vector meson of polarization The mean field approx- 
imation is used here, and the matrix elements needed to evaluate Eqs. (12) and ([H|) can 
be obtained using standard manipulations on Eqs. @ and @. The detailed expressions for 
f s , f v have the same form as those of Ref. |15| in which the nucleus was treated as an ele- 



mentary spherical source of scalar and vector mesons. The difference is that self-consistent 
field theory is used here to obtain the nucleon densities p s ,B- The scalar mesons are found 
to carry 4.70 MeV of the plus-momentum, while the vector mesons carry 253.30 MeV (or 
27%) of the nuclear plus-momentum. The technical reason for the difference is that the 
evaluation of a'(\n, u)a(k, u) counts vector mesons "in the air" and the resulting expression 
contains polarization vectors that give factor of which enhance the distribution of vector 
mesons of low k + . The results for the vector meson distribution are shown in Fig. 2. 

It is worthwhile to see how the present results are related to lept on- nucleus deep inelastic 
scattering experiments. We find that the nucleons carry only 72% of the plus-momentum. 
The use of our in standard convolution formulae lead to a reduction in the nuclear 
structure function that is far too large (~95% is needed |jj) to account for the reduction 
observed in the vicinity of x ~ 0.5. The reason for this is that the quantity M + g s <ft acts 
as a nucleon effective mass of about 670 MeV, which is very small. A similar difficulty occurs 
in the (e, e') reaction JT6| when the mean field theory is used for the initial and final states. 



The use of a small effective mass and a large vector potential enables a simple reproduction 
of the nuclear spin orbit force M. However, effects beyond the mean field may lead to a 



significant effective tensor coupling of the isoscalar vector meson |17[ and to an increased 



value of the effective mass. Such effects are incorporated in Bruckner theory and a light 



front version j0| could be applied to finite nuclei with better success in reproducing the 
data. 

In any case, these kinds of nuclear physics calculations can be done in a manner in which 
modern nuclear dynamics is respected, boost invariance in the z-direction is preserved, and 
in which the rotational invariance so necessary to understanding the basic features of nuclei 
is maintained. 
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TABLES 

TABLE I. Comparison of the single particle spectra of 40 Ca in the equal time (ET) formalism 
(e n — Mtv) with the light front (LF) method (p~/2 — Mw). Scalar and vector meson parameters 
are taken from Horowitz and Serot [9] , and we have assumed isospin symmetry. 



ET 



LF 



State n 

Osi/2 

OP3/2 

Opi/2 



Od 



5/2 



Is 



1/2 



Od 



3/2 



e n - M N (MeV) 
-55.40 
-38.90 
-33.18 
-22.75 
-14.39 
-13.87 



J z = 1/2 

-55.39 

-38.90 

-33.18 

-22.75 

-14.36 

-13.87 



Jz = 3/2 



-38.90 



-22.75 



-13.88 



Jz = 5/2 



-22.74 
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FIGURES 

FIG. 1. Nucleon plus-momentum distribution function, /7v(y), for 40 Ca. Here y = p + /(Ma/A) 



FIG. 2. Vector meson plus-momentum distribution, yf v (y), for 40 Ca. 
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